Centrality Dependence of Multiplicity, Transverse Energy, and Elliptic Flow from 

Hydrodynamics 

P.F. Kolb°'^ U. Heinz", P. Huovinen^ K.J. Eskola'''^ and K. Tuominen'^ 
"■Department of Physics, The Ohio State University, 174 West 18th Avenue, Columbus, OH 43210, USA 
^Institut fiir Theoretische Physik, Universitat Regensburg, D-93040 Regensburg, Germany 
'^Lawrence Berkeley National Laboratory, Berkeley CA 94720, USA 
Department of Physics, University of Jyvaskyld, P.O. Box 35, FIN-40351 Jyvaskyla, Finland 
''Helsinki Institute of Physics, P.O. Box 64, FIN-00014 University of Helsinki, Finland 
(Submitted to Nucl. Phys. A on March 23, 2001) 



O 
O 



m 

cn 
> 

m 
(N 
cn 
o 



Oh' 
I 

Oh: 

CD 



X 



The centrality dependence of the charged multiplicity, 
transverse energy, and elliptic flow coefficient is studied in 
a hydrodynamic model, using a variety of different initializa- 
tions which model the initial energy or entropy production 
process as a hard or soft process, respectively. While the 
charged multiplicity depends strongly on the chosen initial- 
ization, the pT-integrated elliptic flow for charged particles 
as a function of charged particle multiplicity and the pt- 
differential elliptic flow for charged particles in minimum bias 
events turn out to be almost independent of the initial energy 
density profile. 

PACS numbers: 25.75-q, 25.75.Ld 
Keywords: Relativistic heavy-ion collisions; Elliptic 
flow; Hydrodynamic model 

I. INTRODUCTION 

Elliptic flow is a collective flow pattern which de- 
velops in non-central relativistic heavy-ion collisions as 
a result of the spatial deformation of the initial trans- 
verse overlap area. It requires rescattering among the 
produced particles as a mechanism to map the initial 
spatial deformation of the reaction zone onto the finally 
observed hadron momentum distributions. It is quan- 
tified by the second harmonic coefficient V2 of an az- 
imuthal Fourier decomposition of the measured spectrum 
dN / (dy pTclpTdcj)) j^j. Its magnitude V2 and its shape 
V2{pt) as a function of the hadron transverse momentum 
are sensitive to the scattering rate among the produced 
secondaries, especially during the dense early stage of 
the expansion [^j4|. The largest elliptic flow signal, es- 
pecially at high pt, arises in hydrodynamic simulations 
10,^1^ which assume local thermal equilibrium at every 
space-time point, i.e. essentially instantaneous thcrmal- 
ization or infinite scattering rate. Surprisingly, such hy- 
drodynamic simulations |^] are in very good agreement 
with first results from y/s — 130 A GcV Au+Au collisions 
at the Relativistic Heavy Ion CoUider (RHIC) up 
to transverse momenta of 1.5 — 2 GeV/c. 

It was recently suggested P,p|p^ that a combined ana- 
lysis of the full set of hadronic single particle spectra 
and their elliptic flow as a function of collision central- 
ity should allow to outline the domain of applicability 
of the hydrodynamic approach. Inside this domain such 



an analysis would constrain the initial baryon and en- 
ergy density and the flnal freeze-out temperature suffi- 
ciently well to become sensitive to details of the equation 
of state (EOS) of the fireball matter In this context 
V2 provides access to the EOS during the early expansion 
stage which, at RHIC energies, is presumably in the 
quark-gluon plasma (QGP) phase. Direct verification of 
the phase transition between QGP at high energy den- 
sity and a hadron resonance gas at lower energy density 
should then become possible by accurately measuring the 
excitation function of radial and elliptic flow ||]. 

In the hydrodynamic limit the EOS affects the elliptic 
flow signal through the velocity of sound, Cs = \J dP/de 
||l| . However, the sensitivity of V2 on Cg is not very strong, 
and even a first order phase transition, where Cg vanishes 
in the mixed phase, affects V2 only on the 20% level [||. 
This makes accurate measurements and systematic theo- 
retical checks indispensable. One possible source of am- 
biguity, which has not been systematically investigated 
in previous studies [p]J^,|6|,p|,p|,p^ , is the sensitivity of the 
radial and elliptic flow pattern on the shape of the ini- 
tial transverse density proflle. The latter depends on the 
scaling of secondary particle production with the num- 
ber of colliding nucleons which itself is controlled by the 
collision centrality. It was recently found that the 
momentum-space anisotropy at high px, where it is not 
of hydrodynamic origin but due to quark energy loss []l4| , 
is quite sensitive to the initial density profile. A good hy- 
drodynamic baseline for the dependence of f2(pT) on the 
initial profiles may thus help to clarify the physics con- 
trolling the breakdown of the hydrodynamic model and 
the transition to the hard physics domain at high px- 

In the present work, we investigate five options which 
are expected to span the realistic range of possibilities. 
In the first four, either the initial energy or the initial en- 
tropy density are assumed to scale with either the num- 
ber of wounded nucleons ("soft" or "non-perturbative" 
particle production) or the number of binary nucleon- 
nucleon collisions ( "hard" or "perturbative" particle pro- 
duction). In the fifth model perturbative particle pro- 
duction is modified by implementing gluon shadowing 
in the initial state and by limiting the growth of 
the production cross section by gluon saturation in the 
final state ("saturation model" [|6|jl^). This brings in 
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some non-perturbative elements as well. In the first four 
parametrizations we normalize the initial energy density 
profile such that for central collisions we reproduce within 
errors the total charged multiplicity density at midrapid- 
ity, dNc\y/drj\^ri\<i = 555±10%, as measured by PHOBOS 
for Au+Au at V^= ISOAGeV On the other hand, 
in the saturation model the initial energy den- 

sity is fixed by assuming that the transverse energy of 
produced minijets is entirely converted to thermalized 
energy density. The different scaling laws implied by the 
models then translate into different centrality dependen- 
cies of dNch/dr] which can be tested against the new data 
presented recently |l9| , ^ . 

The hydrodynamic model describes an adiabatic evo- 
lution from one local equilibrium state to another. Our 
version of the model assumes longitudinal boost invari- 
ance which implies conservation of the entropy rapidity 
density dS/dy. Using the relation between entropy and 
particle multiplicity in a thermalized system, the mea- 
sured centrality dependence of the final multiplicity den- 
sity dN/dy can thus be mapped onto the centrality de- 
pendence of the initial parton multiplicity density at the 
point of thermalization. This allows to constrain models 
for the initial production of secondary particles, under 
the assumption of subsequent adiabatic evolution. 

We show that different models for the initial energy 
and entropy production lead to different radial and ellip- 
tic flow patterns. At fixed impact parameter, these result 
in different px-dependences of the spectra and of the dif- 
ferential elliptic flow 1)2 (pt)- However, the px-integrated 
elliptic flow as a function of the final charged multi- 
plicity density dNch/dy and the differential elliptic flow 
V2{pt) for minimum bias events show surprisingly little 
sensitivity to the model used for initializing the hydrody- 
namic evolution. These two observables have been used 
to argue for the applicability of the hydrodynamic model 
at RHIC energies [0j8|, and they were shown to be sen- 
sitive to the EOS of the expanding matter HJ^. Their 
insensitivity to the shape of the initial transverse energy 
density proflle may prove to be crucial for the process 
of extracting reliable information on the EOS from a de- 
tailed analysis of the measured collective flow patterns. 

II. MODELS FOR THE INITIAL TRANSVERSE 
ENERGY DENSITY PROFILE 

A. Soft particle production 

In Pb-|-Pb collisions at the SPS, the rapidity densi- 
ties at midrapidity of both the total produced trans- 
verse energy]^ dET/dy, and of the charged multiplicity. 



dNch/dy, scale approximately linearly with the number 
of participating nucleons, A'pait except for very pe- 
ripheral collisions with iVpart~100. Similar earlier obser- 
vations in smaller collision systems and/or at lower en- 
ergies have led to the notion that particle and transverse 
energy production can be described by the "Wounded 
Nucleon Model" |2^] in which each participating nu- 
cleon contributes to the multiplicity and transverse en- 
ergy only once in its first collision and not every time it 
suffers further collisions with other projectile or target 
nucleons. The validity of this model requires destructive 
interference effects between subsequent nucleon-nucleon 
collisions which are thought to be characteristic of non- 
perturbative or "soft" particle production processes. For 
two nuclei colliding at impact parameter 6 = 6 e^. , the 
density of wounded nucleons in the transverse plane 
(parametrized by s = {x,y)) can be calculated from the 
simple geometric formula (Glauber ansatz, see ^3 24 
and references therein): 



7iwN(s;b) = r^(s+i6) 



B 



(1) 



laOOpr; 

1000 
^8 800 
« 600 

Q. 

^ 400- 
200 



part 



10 



12 14 
b{fm) 



FIG. 1. Number of participating ("wounded") nucleons 
and of binary nucleon-nucleon collisions as functions of im- 
pact parameter. This and all following figures refer to Au-|-Au 
collisions at \fs = ISOAGeV. 

Ta is the nuclear thickness function of nucleus A, 

/ + 00 
dzpA{s,z), (2) 
-oo 



with the density pA given by a Woods-Saxon profile. 



PA{r) 



1 + cxp[(r - i?o)/e] 



(3) 



^We use Et ~ Ei where the sum is over all particles. 



and similarly for nucleus B. For the Woods-Saxon pro- 
file we take standard parameter values, Rq = 1.12 A^/'^ — 
0.86 A~^/'^fm for the radius and ^ = 0.54fm for the sur- 
face diffuseness. The integral over (|^) gives the number 
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of "wounded" or "participating" nucleons, A^'part = -/Vwn, 
as a function of impact parameter (see Fig. p. For the 
total nucleon-nucleon cross section at -^s — 130 A GeV we 
take a = 40 mb. 



1. Model sWN 

Hydrodynamics with boost-invariant longitudinal ex- 
pansion conserves the entropy per unit rapidity dS/dy. 
At fixed freeze-out temperature and chemical poten- 
tial, dS/dy is related one-to-one to the measured 
charged multiplicity density, dN^h/dy. If the latter 
scales linearly with the number of wounded nucleons, 
A^wn(&) = / rf^s ^wn(s; b), barring a strong 6-depen- 
dence of the freeze-out conditions it is natural to assume 
that the initial entropy density in the transverse plane is 
proportional to that of the wounded nucleons [l|j25[|: 



; To; b) = KsItq) nwN(s; b) . 



(4) 



At SPS and RHIC energies the net baryon density is so 
small that its influence on the pressure and thus on the 
developing flow pattern is hardly noticeable. For simpli- 
city we can thus also parametrize the initial net baryon 
density n as being proportional to nwN: 



n{s; To; b) = X„(ro) nwN(s; b) . 



(5) 



Equations and (|^) together define the first model for 
the initialization, sWN, to be studied below. 

The implementation of the initialization into the 
hydrodynamic computer code is complicated by the fact 
that in our formulation the hydrodynamic equations 
propagate the fields n{x) and e(x) (and not s{x)). The 
initialization (^) thus requires the additional step of com- 
puting e{s]To]b) from s{s\TQ]b) by solving an implicit 
equation resulting from the thermodynamic identity 



r(e, n) s — e + p(e, n) — /i(e, n) n , 



(6) 



where T(e, n), p(e, n), and /i(e, n) are tabulated values for 
the temperature, pressure and baryon chemical potential 
corresponding to the selected equation of state. 



quickly with A^wN than dNch/dy. There is some indica- 
tion for this to happen in the WA98 data from Pb-f Pb 
collisions at the SPS where E^/Nch at midrapidity 
rises slightly from very peripheral to semiperipheral col- 
lisions and then saturates from semiperipheral to central 
collisions. However, the experimental effect is small and 
not clearly statistically significant. 

Given this unclear experimental situation, and because 
it simplifies the initialization process, we assumed in our 
previous work that it is the initial energy den- 

sity (and not the entropy density) which scales with the 
density of wounded nucleons in the transverse plane: 



e(s; To; b) = Ke{To) n-^^{s; b) . 



(7) 



Equation (0) together with Eq. (||) defines our second 
initialization model, eWN. The effect of flow on E^/Nch 
at RHIC and its dependence on the initialization of the 
hydrodynamic evolution will be discussed in Section III. 



B. Hard particle production 

At higher and higher collision energies, one expects 
that hard collisions among quarks and gluons from the 
nuclear structure functions become more and more im- 
portant and eventually dominate secondary particle pro- 
duction p^ ]. In this limit secondary particle produc- 
tion is a result of incoherent parton-parton collisions, and 
each nucleon-nucleon collision contributes equally to the 
cross section. Particle and energy production should then 
be related to the distribution of the number of binary 
nucleon-nucleon collisions in the transverse plane. It is 
given in terms of the nuclear thickness functions (|^) by 



dNt^{s;b) 



= aTA{s+lb)TB{s-\b) , 



(8) 



where a is the nucleon-nucleon cross section. The in- 
tegral of (H) is the nuclear overlap function, TAB{b) = 
J d^sTA (s-t-ifo) Tb (s—^b), which is normahzed to AB: 
J d'^hTAsib) = AB. It gives the total number of binary 
collisions as a function of impact parameter, Nco\\{b)^ 
which is also shown in Fig. |l]. 



2. Model eWN 



1. Model sBC 



In a hydrodynamic approach with fixed freeze-out con- 
ditions, it is strictly speaking not possible that both the 
transverse energy dE^ / dy and the charged multiplicity 
dNch/dy scale linearly with A^wn- During the hydrody- 
namic evolution thermal energy is converted into longitu- 
dinal and transverse flow energy, and the time available 
for this conversion increases with the number of wounded 
nucleons. Since the transverse flow increases the average 
transverse energy per particle, dE^/dy should rise more 



If the system of secondary particles thermalizes quickly 
via elastic collisions, their number density defines, up to 
a proportionality constant, the initial entropy density at 
the beginning of the hydrodynamic expansion: 



s[s;To;b) = ks{To)TA {s+^b)TB {s- 



(9) 



Assuming that the initial net baryon density in the trans- 
verse plane can also be calculated perturbatively from the 
nuclear structure functions ETI , we write 
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n(s; To; 6) = 7^„(ro) Ta {s+\b) Tb 



(10) 



As discussed in Section |II A l| , the results do not de- 
pend on whether we use ([lOD or (H) . The combination 
of Eqs. and ( p^ defines model sBC for the initializa- 
tion. 



2. Model eBC 

Within the perturbative approach to particle produc- 
tion one can also argue that each nucleon-nucleon colli- 
sion contributes with equal probability not only to the 
number of produced secondaries, but also to the energy 
carried by them. This leads to the ansatz that the initial 
energy (and not the entropy) density is proportional to 
the density of binary collisions in the transverse plane: 



e(s; To; b) = K,{to) Ta (s+^b) Tb (s-^b) 



(11) 



This equation together with Eq. ( |l0| ) defines model eBC 
for the initialization. 



C. The saturation model 

In the saturation model [ p^p7| ] the production of par- 
tons becomes inhibited below a saturation scale Psat, de- 
termined as the transverse momentum scale where the 
produced partons start to overlap transversally. The for- 
mation time of the QGP at each point s is given by 
Tsat = l/Psat(-\/s, A, B, s; b) . The local initial energy den- 
sity profile can then be worked out at the central slice as 
in I2II, 



e{s-b) 



cPs dz 
TA{s+\b)TB{s-\b) 

xa{ET){^/^,Pse,t,Ay,A,B,s;b) 



(12) 



where dz « TgatAy within the central rapidity unit Ay. 
The first i^x-moment of the minijet distribution, a{ET) 
p6| , p7| , is computed in lowest order pQCD as the first pt- 
moment of the distributions of minijets with > Psat 
and y within Ay. The EKS98 shadowing |l^] of the par- 
ton distributions is included and the NLO contributions 
to (j{Et^) pofl are simulated by an approximate factor 
if = 2, as also done in [|l6|,0. 

In |l^ a fully saturated system was considered by 
extending the computation down to very low values of 
Psat = 0.5 GeV, and neglecting the tails of the number 
(energy) density distributions at large transverse dis- 
tances. Now, for the hydrodynamic description, we have 
to consider the tails also. To maintain the spirit of the 
local saturation model of ||l^, we restrict the satura- 
tion to the regions where psat > 0.75 GeV and compute 



the tail of the energy density profile from Eq. (|lj) with 
Psat = 0.75 GeV. Through this procedure, the multiplici- 
ty of 1 1^ in central collisions is recovered (if computed as 
in ||l7|). For central collisions, the tail contributes only 
7% to the total multiplicity. We would like to emphasize 
that we have not tried to fit the saturation model to the 
centrality data but, rather, to keep the approach as close 
to the orginal local saturation idea as possible. 

As discussed in [|l6| , from the energy per particle point 
of view the system looks thermal already at saturation. 
The same can be shown to hold also in the locally sat- 
urated system |l^. Within the saturation framework it 
is therefore meaningful to switch on the hydrodynamic 
evolution already at Tgat which now is a local quantity. 
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FIG. 2. The formation and initial thermalization time 
Tsat (&) = Tsat (s = 0; 6) = 1/psat in the saturation model, as a 
function of impact parameter 6. 

For the hydrodynamic set-up used in this study, how- 
ever, the density profiles at a constant initial time are 
required. To circumvent this problem, we consider two 
possibilities: first, for each impact parameter b, we sim- 
ply use the earliest time, i.e. t{s = 0), as the initial time. 
These are shown in Fig. ^j. Notice that the maximum ini- 
tial time is (0.75 GeV)"^ — 0.26 fm/c, and that a constant 
initial time is a good approximation for the central re- 
gion and for small impact parameters. Alternatively, the 
energy densities computed from Eq. ( p^ ) can be evolved 
at each s and at each Tsat(s;6) to Ti = 0.26 fm/c assum- 
ing boost-invariant flow without transverse expansion:^ 
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We 



e{s;n;b) = e(s; Tsat(s; b); b) x (rsat(s; b)/Ti) 
have checked, however, that the latter procedure leads 
only to a few percent increase in the central multiplici- 
ties relative to the former one. Therefore, we take in 
the following the initial time for the saturation model 
from Fig. and use e(s; Tsat(&); b) as the initial profile 



■^We thank P.V. Ruuskanen for helpful discussions related to 
this point. 
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for the transversally expanding hydrodynamics, as dis- 
cussed next. 



D. Initial energy density and spatial anisotropy 

Figure |3| shows the initial energy density profiles for 
the different initialization models, {x is the direction 
of the impact parameter b, while y points perpendicu- 
lar to the reaction plane.) For models sWN, eWN, sBC 
and eBC the profiles are normalized such that in each 
case the total charged multiplicity per unit pseudorapid- 
ity at 6 = is dN^Y^/dy = 670 at y = 0. For model eWN |] 
this corresponds to (iiVch/d?7||,,|<i =545 for the 6% most 
central collisions (a bit less for the three other models). 
Within errors, this is consistent with the first published 
PHOBOS data [Q, but slightly below the more accurate 
recent data from PHENIX and PHOBOS [|l9|,|o| . To pre- 
serve consistency with our previous publications [^j^ , we 
decided against increasing our normalization of dNch/drj 
to the new data since this would have implied retuning 
the initial conditions and freeze-out temperature in order 
to keep the spectra and elliptic flow unchanged. 
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FIG. 3. Initial energy density profiles in the transverse 
plane for the five different initialization models described in 
this Section. Top: 6 = 0. Bottom: 6 = 10fm. In each case 
two cuts in a;-direction are shown, one for 2/ = and one for 
i/ = 5fm. For the saturation model ("sat.") the profile was 
hydrodynamically propagated from the initial thermalization 
time assumed in that model to ro = 0.6 fm/c where the hy- 
drodynamic evolution was started for the other four models. 



For the comparison with the other models in Fig. ||, 
the profile obtained from the saturation model was 
evolved from the formation time Tsat(&) (see Fig. ^) to 
tq = 0.6 fm/c, the time at which the hydrodynamic evolu- 
tion was started for the other four initializations. Within 
the line widths of Fig. ^, it did not matter whether this 
scaling was done by assuming only boost-invariant lon- 
gitudinal expansion, e(s,To) =e(s,Tsat) x (rgat/To)'*/^, or 
by including also the transverse expansion. 

Since for a thermalized parton gas at nearly vanish- 
ing net baryon density e~s'*^'^, Eq. (^) gives distribu- 
tions which are more sharply peaked at the origin than 
Eq. (prl). We will see that this results in a differ- 
ent centrality dependence of the total entropy per unit 
rapidity dSjdy (which, for boost-invariant longitudi- 
nal expansion, is proportional to the integral of s(s) 
over the transverse plane). Model eBC interpolates be- 
tween model sWN, where dS/dy ^ Np^-^t, and model sBC, 



where dS/dy ^ NcoW ^ -^part (the latter proportionality 
was checked numerically to hold with excellent accuracy 
over the entire impact parameter range). This is similar 
to the HIJING model |^ where the final charged par- 
ticle rapidity density is a linear superposition of a soft 
component, scaling with A'part, and a hard component 



which scales with Nc, 
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In contrast to the mo- 



VcoU ^'part- 

dels here, however, HIJING has no rescattering among 
the produced particles and thus no collective flow. Fi- 
nally, the saturation model is seen to be close to model 
eWN for central collisions, while for semiperipheral col- 
lisions it nearly coincides with model sWN. 
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FIG. 4. Initial spatial anisotropy as a function of impact 
parameter, for the different initializations. 



Figure |j shows the initial spatial anisotropy 



e4b) = Sib) 



(13) 



as a function of impact parameter, evaluated with the ini- 
tial transverse energy density as weight function, for the 
five initialization models studied here. At fixed impact 
parameter, varies by up to 20%. Since the hydrody- 
namic evolution maps the initial spatial anisotropy onto 
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the final momentum-space, one expects from Fig. ^ that 
the impact parameter dependence of the elliptic flow V2 
should show a similar sensitivity to the initialization. We 
will see that, except for the saturation model, this model 
sensitivity is almost exactly cancelled by the correspond- 
ing variations in the impact parameter dependence of the 
produced charged particle multiplicity (i.e. entropy). 



III. CENTRALITY DEPENDENCE OF 
MULTIPLICITY AND TRANSVERSE ENERGY 

Using the initial transverse energy and baryon density 
profiles from the previous section and the numerical code 
described in |P , we solve the hydrodynamic equations for 
the transverse evolution of the reaction zone, assuming 
boost-invariant expansion in the longitudinal direction. 
We study only Au+Au coUisions at ^/s = 130 A GeV. For 
the models eWN, sWN, eBC and sBC wc use the same 
initial time tq = 0.6fm/c as in to start the hydrody- 
namic expansion. For the saturation model the hydrody- 
namic expansion is started at Tgat (b) shown in Fig. ^. All 
calculations reported here are performed with EOS Q, an 
equation of state with a first order phase transition from 
a hadron resonance gas to a non-interacting quark-gluon 
gas at Tc = 164 MeV (Sl), and freeze-out at Tf = 120 McV. 
The sensitivity of the spectra and elliptic flow to the EOS 
and freeze-out temperature were studied in Mu ■ 



A. Charged particle multiplicity 

Figure |^ shows the final charged particle rapidity den- 
sity per participating nucleon pair resulting from the dif- 
ferent initializations. In the bottom panel the rapidi- 
ty density at midrapidity is converted to pseudorapidity 
density by the transformation 



drj 



E 

z£ charged 



Pt 



(14) 



setting y=pi^—0. At fixed freeze-out temperature 
dNch/dy is a direct measure for the entropy density 
dS/dy, and the approximate constancy of the curve for 
model sWN thus reflects the approximate conservation 
of entropy by the hydrodynamic evolution. In fact, 
the slight increase of {dNch / dy) / Np^it for small values 
of iVpart (Fig. H, upper panel) can be traced back to 
a small amount of entropy production by deflagration 
shocks which arise during the hydrodynamic expansion 
stage as a consequence of the first order phase transition 
p2[ once the initial energy density in the center of the 
reaction zone increases above the critical value for QGP 
formation. 

Model eWN (solid line in Fig. ||) was used in Refs. 
|,||J|, while model sWN was employed in |l|,|l2|,||. 



Fig. H shows that both initialization models are dis- 
favoured by the data, model eWN more so than mo- 
del sWN. The saturation model produces a rather simi- 
lar centrality dependence of the charged multiplicity as 
model eWN. 
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FIG. 5. Charged particle yield per participating nucleon 
pair at midrapidity as a function of the number of partici- 
pants. All curves were normalized to dN^h/dy = 670 for cen- 
tral (6 = 0) collisions (see discussion below Fig. ^). The top 
panel gives the rapidity, the bottom panel the pseudorapidity 
density. The data are taken from Refs. jl^-|2o|]. 

As expected, the saturation model results are close to 
those obtained in Ref. |l^. The careful reader notices, 
however, that the values ((iA'ch/(iy)/(0.5A^part) for the 
saturation model shown in Fig. || are 0(10%) lower than 
those shown in Fig. 4 of [OL This is due to a number of 
partially cancelling effectspSf : (i) Inclusion of the quark 
degrees of freedom increases the initial entropy by about 
20% compared to |l^ . (ii) In a realistic hadron resonance 
gas containing also heavy particles, at a freeze-out tem- 
perature of T[ = 120 MeV each particle carries on average 
4.87 units of entropy, instead of the 4 units assumed in 
JlTf (which is a good approximation in the ultrarclativis- 
tic limit) . This decreases the total multiplicity from that 
in (iii) Finally, the decay of unstable resonances 

changes the fraction of charged particles from 2/3 to 
w 0.62 of the total multiplicity. The net reduction effect 
for the final charged multiplicity rapidity density is about 
10% as mentioned. 
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The best representation of the shape of the data, even 
down to rather peripheral collisions (low TVpait-values), 
seems to be given by model eBC, once appropriately 
renormalized to the new, slightly higher multiplicity den- 
sities in central collisions. Alternatively, one may try a 
linear combination of the sWM and sBC parametriza- 
tions, with a large sWN and a smaller sBC contribu- 
tion (the latter is apparently absent at the SPS 
where the multiplicity scales linearly with A'part)- As 
discussed in Section II C , the main difference between 



the model eBC and the saturation model is that the 
s-dependence in eBC is contained only in the product 
Ta {s+^b^ Tb (s— ifo) whereas in the saturation model 
[0 also the cross sections (and their moments) depend 
on s through the local saturation scale Psat(^)- The data 
do not appear to support this particular implementation 
[p7| of gluon saturation at the present RHIC energy of 
y/s = 130 A GeV for the non-central collisions. In central 
collisions, however, parton saturation is not excluded as 
the dominant mechanism. We conclude that obviously 
the saturation region is now extended too far in trans- 
verse direction and that the eBC tails in the initial energy 



density profile (which were not included at all in |17|) 
should be given more weight than they have now. On the 
other hand, the data strongly indicate that the initial en- 
ergy deposition process does involve a component which 
scales with the number of binary nucleon-nucleon colli- 
sions. We refrain from a discussion whether or not this 
indeed proves the onset of "hard" perturbative physics 
at RHIC ioi. 



B. Transverse energy per particle 

Hydrodynamic flow is a result of the conversion of a 
fraction of the thermal energy into collective flow ki- 
netic energy, through work done by the thermodynamic 
pressure. As a consequence, the thermal energy per 
particle decreases as a function of time. If the reac- 
tion zone were to undergo boost-invariant expansion only 
along the longitudinal direction, dS/dy ~ dNch/dy would 
be constant and dE^/dy would be proportional to the 
thermal energy density. The reduced thermal energy 
per particle would thus be directly reflected in the fi- 
nally observed transverse energy per charged particle, 
(dET/dy)/ (dNch/dy) Since in more central colli- 

sions the thermodynamic pressure does longitudinal hy- 
drodynamic work for a longer time, the finally observed 
ratio (dE^/dy)/ (dNch/dy) should increase much more 
slowly with A'^part than the initial energy per particle | |3^ ] 
established during the energy deposition process. In real- 
ity, the system undergoes not only longitudinal, but also 
transverse expansion. The transverse collective flow adds 
a kinetic contribution to E^ which reduces the diluting 
effect on E^/Nch from the longitudinal expansion. 



This shows that the centrality dependence of E^/Nch 
is an interesting observable which reflects the interplay 
of the centrality dependences of initial transverse energy 
and particle production, longitudinal work done by the 
hydrodynamic pressure, and radial flow. In Figure |^ we 
show the corresponding results from the hydrodynamic 
model for Au-hAu collisions at RHIC (Vs = 13G^GeV) 
for the five different initialization models. For compari- 
son we also show the centrality dependence for Pb-|-Pb 
collisions at the SPS, for model eWN with parameters 
adjusted to the measured spectra |5||^. Comparing the 
eWN curves for SPS and RHIC, one sees that the larger 
radial fiow at RHIC has very little effect on the average 
i?T per particle in the final state: the gain in transverse 
kinetic flow energy is largely compensated by more longi- 
tudinal work done at RHIC (smaller tq , somewhat larger 
Tf). The eWN curve at the SPS is consistent in shape 
with the WA98 data shown in Fig. 14 of but at the 
upper edge of their systematic error band. 
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FIG. 6. Average transverse energy per charged particle at 
midrapidity as a function of the number of participants. Re- 
sults are shown for all five initialization models for Au-f-Au 
collisions at y/s = 130 A GeV. In addition we plot the result 
from model eWN for Pb-|-Pb collisions at = 17 A GeV, us- 
ing the initial conditions from Ref. [^,^ (thin dashed line). 

At RHIC, the models eWN and sWN as well as the 
saturation model all give similarly flat centrality depen- 
dences for the average Et per charged particle; they 
differ slightly in normalization and in how steeply this 
ratio drops for very peripheral collisions. The data in 
Figure ^ disfavour these initialization models. On the 
other hand. Model eBC, which works better in Figure ||, 
exhibits a markedly stronger centrality dependence of 
E^/Nch, with a clearly visible slope below A^part — 250. 
(This is even more true for the initialization sBC, but 
for that model the corresponding curve in Fig. || is too 
steep.) This reflects a similar strong centrality depen- 
dence of [dE^ / dy) / (dS / dy) in the initial state, which is 
only partially offset by the increasing amount of longitu- 
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dinal hydrodynamic work done as the collisions become 
more central. 

As recently pointed out in | |3^ , a measurement of the 
ratio plotted in Fig. ^ will be very helpful in comple- 
menting the information contained in Fig. || and check- 
ing for the effects of longitudinal work. We note that 
non-equilibrium expansion scenarios like those studied in 
p5[ , where the matter performs less longitudinal work 
than predicted hydrodynamically, should lead to an even 
stronger dependence of /Nch on iVpart than that shown 
in Fig. |6[ In contrast, recent data from the PHENIX Col- 
laboration |36j exhibit the opposite tendency: they are 
flatter than the eBC curve in Fig. |^ and, in fact, quantita- 
tively consistent in magnitude and shape with the WA98 
data at the SPS Q. It remains to be seen whether e.g. 
a linear combination of a large sWN with a small sBC 
component (which can be made to agree with the mul- 
tiphcity data in Fig. || Q) leads to a sufHciently flat 
behaviour in Fig. || to be consistent with the PHENIX 
data 1 36 1, and whether the magnitude of E^/Nch can be 
lowered by replacing the unrealistic chemical composition 
from the hydrodynamical model by a more realistic one 
(our hydrodynamic model assumes that chemical equilib- 
rium is preserved until kinetic freeze-out at Tf = 120 MeV 
whereas preliminary STAR data indicate a much higher 
chemical freeze-out temperature of around 170-180 MeV 
0)- 



IV. RADIAL AND ELLIPTIC FLOW 
A. Radial flow 

The single particle spectra for positive pions and an- 
tiprotons for the different initialization models are shown 
in Figure 0, for two representative values of the impact 
parameter. For central (6 = 0) collisions, the systematics 
of the mx-slopes shows that for models eWN, sWN, eBC 
and sBC the amount of radial flow created in the collision 
increases in the same order as the central initial energy 
density shown in Figure |^. The saturation model is an 
exception: it gives larger radial flow than both eWN and 
sWN although, at tq = 0.6 fm/c, its central energy density 
is even below that of model eWN. The reason is that for 
the saturation model the hydrodynamic evolution starts 
earlier and the build-up of transverse flow thus begins 
already at time Tsat(^) (see Fig. ||). 

The effect of the initialization on the single particle 
spectra is significant, and switching from model eWN 
used in |^ to model eBC, which according to Figure ^ 
provides a better fit to the centrality dependence of the 
charged multiplicity, considerably flattens the single par- 
ticle spectra predicted in Ref. |^ for central collisions. 
This model sensitivity is weaker for semiperipheral colli- 
sions (6 = 8fm), again in agreement with Figure || which 
shows that in this range the initial energy density profiles 



are rather similar for all the models. If hydrodynamics 
still applies to even more peripheral collisions, models 
eBC and sBC would there predict steeper spectra than 
models eWN and sWN. 
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FIG. 7. Transverse mass spectra for positive pions and an- 
tiprotons, for central (upper set of curves) and semiperipheral 
collisions (lower set of curves, divided by 10 for clarity), for 
the five initialization models studied in this paper. 



B. Elliptic flow 

Figure ^ shows that for the four models eWN, sWN, 
eBC, and sBC the differential elliptic flow V2{pt) fol- 
lows a similar pattern as the radial flow: the model with 
the largest initial central energy density produces the 
strongest elliptic flow. The saturation model docs not 
follow this systematics, which again is presumably due 




2 2 S 

Pt (GeV) 

FIG. 8. The px- differential elliptic flow V2(pt) for charged 
hadrons, at three representative impact parameters and for 
the five initialization models studied in this paper. 

to the earlier start of the hydrodynamic evolution in 
this model. The impact parameter dependence of this 
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model sensitivity is now controlled by the initial spa- 
tial anisotropy ex- According to Figure ^, for most of 
the impact parameter range it satisfies the hierarchy 
e,(eWN) < e,(sWN) < e,(sat) < e,(eBC) < e,(sBC), 
and this is clearly reflected in the slopes of V2{pt) for 
Pt < 2GeV/c, as shown in Figure ||. Again, the effect 
is significant and can reach up to 20%. At large impact 
parameters, b > 10 fm, and for pT > 2GeV/c this simple 
ordering is broken. 

In view of these results, it is somewhat of a pleasant 
surprise to see in Figure ^ that the px-integrated ellip- 
tic flow V2 for all charged hadrons, when plotted as a 
function of the midrapidity density of charged particle 
multiplicity as an experimentally easily accessible cen- 
trality measure, shows very little sensitivity to the initial- 
ization model. The models eWN, sWN, eBC, and sBC 
yield almost identical results, and only for the saturation 
model V2 is slightly larger (except for the most periphe- 
ral collisions). Comparing the two extreme models eWN 
and sBC, one sees from Figure |] that, at fixed charge 
multiplicity, the latter favours smaller impact parame- 
ters which correspond to smaller and U2(pt) (Figs. ^ 
and but flatter single particle spectra which give more 
weight to larger V2 at higher px (Fig- §|) ■ The net effect of 
this intricate interplay is an almost complete cancellation 
of these counteracting tendencies. 
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FIG. 9. The pT-integrated elliptic flow 112 for charged 
hadrons as a function of midrapidity charged multiplicity den- 
sity, for different initialization models. The data are taken 
from Ref. . (For a discussion of the horizontal axis rich /wmax 
in theory and experiment, we refer to Refs. (0,^.) 

A similar cancellation happens for the differential 
anisotropy W2(pt) for charged hadrons in minimum bias 
events (left panel of Figure |l0|) : 



W2(Pt) 



^ bdb V2{pt; b) 



dy Pt dp-y 



Jbdb 



dNa] 



dy Pt dp^ 



:ib) 



(15) 



At fixed b, the models with smaller elliptic flow are 
weighted with larger charged multiplicities, but slightly 



steeper spectra, and again the net result is an almost 
miraculous cancellation of all sensitivities to the initial- 
ization model for px ^ 1.5 GeV/c. Above px = 1 GeV/c, 
the saturation model gives the largest V2{pt)- 

This cancellation does not carry over to heavier parti- 
cles. The right panel in Figure |l^ shows the differential 
elliptic flow V2{pt) for identifled protons and antiprotons, 
which is seen to exhibit more signiflcant variations as 
the initialization of the hydrodynamic model is changed. 
Generically, the model which gives the largest radial flow 
at small impact parameters produces the smallest proton 
elliptic flow at small values of px. This is qualitatively 
consistent with the general analytic discussion of radial 
flow effects on the elliptic flow for heavy particles pre- 
sented in Ref. §. Again, the saturation model gives the 
largest elliptic fiow of all studied initializations for pro- 
tons with Px > 1 GeV/c. 




Pt (GeV) 
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FIG. 10. The differential elliptic flow V2{pt) for minimum 
bias events, for all charged hadrons (left) and for identified 
protons or antiprotons (right). The curves correspond to dif- 
ferent initialization models as indicated in the Figure. The 
data in the left panel are from Ref. |7| . 



V. CONCLUSIONS 

Within a hydrodynamic model with exact longitudi- 
nal boost invariance, we have shown that the centrality 
dependence of the production of particles, transverse en- 
ergy and transverse flow is significantly infiuenced by the 
shape of the initial energy density profile in the transverse 
plane. This profile is intimately related to the nature of 
the primary particle production process which converts 
beam energy into the matter forming the collision fireball. 
The available data for the charged multiplicity in Au+Au 
collisions at RHIC as a function of collision centrality 
are compatible with an initial energy deposition process 
involving a component which scales with the number 
of binary nucleon-nucleon collisions. They strongly dis- 
favour the wounded nucleon parametrization used in ear- 
lier hydrodynamic simulations |l]j|,^,||,| 
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In non- 
central collisions they are, at the present RHIC energy. 
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also at variance with predictions from the gluon satura- 
tion model ll^. The best agreement with the multiplic- 
ity data at y/s = 130AGeV [|l|-|2^ would be obtained 
by assuming that most of the initial thermalized energy 
density is proportional to the number of binary nucleon- 
nucleon collisions, with a smaller "soft" contribution pro- 
portional to the number of wounded nucleons. If one in- 
stead parametrizes the initial entropy density as a super- 
position of "hard" and "soft" components, as implied by 
the approaches in | ]30| , p3[ , one needs a larger "soft" and 
smaller (~ 10%) "hard" component |33|. A model with 
similar properties is obtainable also in the saturation ap- 
proach considered here by limiting the saturation region 
to smaller transverse distances and thereby increasing the 
contribution from the tail which scales with the number 
of binary collisions. All these models predict an average 
transverse energy per particle in the final state which 
rises slowly with increasing number of participating nu- 
cleons. A measurement of E^/Nch as a function of A^part 
is sensitive to the longitudinal work done during the hy- 
drodynamic evolution and can be used to check the con- 
sistency of the hydrodynamic approach. 

The centrality dependence of the px-averaged elliptic 
flow V2 and the differential elliptic flow 1^2 (pt) for all 
charged particles from minimum bias collision events ex- 
hibit only minor sensitivity to the shape of the initial 
energy density profile. Thus, the conclusions extracted 
from earlier studies of elliptic flow using the wounded 
nucleon parametrization (which is now disfavoured by 
the new multiplicity data) remain valid. For identified 
heavy hadrons such as (anti)protons the differential el- 
liptic flow V2{pt) at low pt is smaller for the binary col- 
lision induced initial conditions than for the wounded nu- 
cleon parametrizations. This stems from the somewhat 
stronger radial flow created by the "hard" initializations, 
which should also be visible in the single-particle spectra. 

In closing we note that the initialization models stud- 
ied in the present paper assume a one-to-one correspon- 
dence between impact parameter, number of participat- 
ing nucleons, and initial energy density proflle in the 
transverse plane. This neglects the possibility of strong 
event-by-event fluctuations in the initial density profile 
at fixed number of charged particles in the final state. 
As pointed out recently p7| , strong event-by-event fluc- 
tuations in the initial energy density profile may result 
in sizeable event-by-event fiuctuations of V2 at fixed im- 
pact parameter. Our present study shows that small and 
smooth variations in the initial energy density profile at 
fixed impact parameter cause effects on the particle mul- 
tiplicity and on the radial and elliptic fiow which compen- 
sate each other in such a way that W2(?^ch/'^max) is almost 
unaffected. It would be interesting to study whether this 
compensation also survives the much larger fluctuations 
studied in Ref. fs^. 
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